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Abstract 

We establish spectral estimates at a critical energy level for /i-pseudors . Via 
a trace formula, we compute the contribution of isolated (non-extremum) 
critical points under a condition of ’’real principal type”. The main result 
holds for all dimensions, for a singularity of any finite order and can be in- 
variantly expressed in term of the geometry of the singularity. When the 
singularities are not integrable on the energy surface the results are signi¬ 
ficative since the order w.r.t. h of the spectral distributions are bigger than 
in the regular setting. 
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1 Introduction. 

If Ph is a self-adjoint /i-pseudo-differential operator, or more generally /^-admissible 
(see El), acting on a dense subset of L 2 , a classical and accessible problem is to 
study the asymptotic behavior, as h tends to 0, of the spectral distributions : 

7 (£>,*.)= £ v(VErT). (1) 

|Aj(/o--e|<s 

where the A j(h) are the eigenvalues of Ph- Here we suppose that the spectrum 
is discrete in [E — e, E + e], a sufficient condition for this is given below. One 
motivation is that it is in general not possible to compute the spectrum and one 
has to use statistical methods to gain spectral information. A second motivation is 
the existence of a duality between spectrum of quantum operators and the classical 
mechanic attached to their symbols. 

If po is the principal symbol of Ph, an energy E is regular when Vpo(x , £) ^ 0 
on the energy surface E e = {(£,£) € T*W n / po(x,£) = E} and critical when it is 
not regular. It is well known that asymptotics of 0 , as h tends to 0, are closely 
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related to the closed trajectories of the classical flow of po on the surface Ev¬ 
idence, there is a duality between the following objects : 

limy (E,ip,h) ^ {(t,x,£) G supp(^) x E s / $t(x,£) = (x,£)}, 

h—> 0 

where <f> t is the flow of the Hamiltonian vector field H po = d^po-d x — d x po-d f. 
This duality has a universal character and does not systematically requires the 
presence of an asymptotic parameter, as can show the trace formulae of Selberg and 
Duistermaat-Guillemin [5]. In the semi-classical setting, this relation was initially 
pointed out in the physic literature DD&D- From a mathematical point of view, 
and E a regular energy, a non-exhaustive list of references concerning this subject 
is 0], 021, TT and more recently with a different approach [7j. See also jl"2| for 
the case of elliptic operators. 

If E is no more a regular value, the behavior of 0 depends on the singularities 
of p on Tie which leads to technical complications. The case of non-degenerate criti¬ 
cal energies, that is such that the critical-set £(po) = {(x,£) G T*W n / dp 0 (x,£) = 
0} is a compact C°° manifold with a Hessian d 2 po transversely non-degenerate 
along this manifold has been studied first by Brummelhuis et al. [J. The problem 
was solved for quite general operators but for some ’’small times”, i.e. for supp(</3) 
contained in a neighborhood of the origin so that the only period of the linearized 
flow in supp((/3) was 0. Later, Khuat-Duy |1SI has obtained the contributions of the 
non-zero periods of the linearized flow with supp(<£) compact, but for Schrodinger 
operators with symbol £ 2 + V ( x ) and a non-degenerate potential V. Our contri¬ 
bution was to generalize the result of m for more general operators but under 
extra assumptions on the flow (see 0). Finally, in [Bj the case of totally degener¬ 
ate extremum was treated and the objective of this work is to study degenerate 
singularities which are not an extremum of the symbol. 

After a reformulation, based on the theory of Fourier integral operators, the 
asymptotics of 0 can be expressed in terms of oscillatory integrals whose phases 
are related to the flow of po on E e- When (xo,£o) is a critical point of po, it is 
well known that the relation Ker(d Ci £$ t (xo, Co) — Id) ^ {0} leads to the study 
of degenerate oscillatory integrals. In this work we consider the case of a totally 
degenerate energy, that is such that the Hessian matrix at our critical point is zero. 
Hence, the linearized flow for such a critical point satisfies d®,£$t(xo, £o) == Id, for 
all t G R. A fortiori : 


Ker(d x ^ t (x 0 ,Zo) - Id) = T X0 ^ 0 T* K" ~ R 2n , Vi G R, (2) 

and the oscillatory integrals we have to consider are totally degenerate. In particu¬ 
lar, it is impossible to use the stationary phase method to determine the asymptotic 
behavior of Eq. 0 . 

To solve this problem we will establish suitable normal forms, for the phase 
functions of our oscillatory integrals, for whom it is possible to generalize the 
stationary phase formula. The construction is geometric and is independent of 
the dimension but the asymptotic expansion of the related oscillatory integrals 
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depends on the dimension and on the order of the singularity at the critical point. 
Finally, since the normal forms have a geometrical meaning, it is possible to express 
invariantly the top order coefficients of the asymptotic expansions in term of the 
geometry of the singularity on the energy surface. 


2 Hypotheses and main result. 

Let Ph = Op™(p(x, £, h)) an /i-pseudodifferential operator, obtained by Weyl quan¬ 
tization, in the class of /i-admissible operators with symbol p(x, £, h) E h j Pj{x,C), 
i.e. there exists pj £ E^T^R") and such that : 

P h =J2 hj Pj ( x ’ hD *) + h N R N (h), VN £ N. 

j<N 

Here is a bounded family of operators on L 2 (R n ) for h < ho and : 

E m( T * R n) = | a . T * R n _ + Cj gup | d^d^a(x,^)\ < C atf} m(x,£), Va,/? £ N”}, 

where to is a tempered weight on T*R". For a detailed exposition on /i-admissible 
operators we refer to the book of Robert ng In particular, po(x, £) is the principal 
symbol of Ph and pi(x, £) the sub-principal symbol. We note d> f := exp (tH Po ), the 
Hamiltonian flow of H po = d^po-d x — d x po-d$ and E^ = p(f 1 (E) C T* R™ the 
energy surfaces of po . 

We study here asymptotics of the spectral distributions : 

7 (E c ,<p,h) = y( Aj ^ E ° )i h-+ 0 + , (3) 

A j (h)^[E c —£,E c -\-e] 

under the hypotheses (Hi) to (H 4 ) given below. We use here the notation E c to 
recall that this energy will be chosen critical. 

(Hi) The symbol of Ph is real and there exists £q > 0 such that the set Pq ([E c — 
£0 ,E C + £ 0 ]) is compact in T* R”. 

Remark 1 By Theorem 3.13 of m the spectrum cr(Ph) H [E c — £,E C + e] is 
discrete and consists in a sequence \i(h) < \- 2 (h) < ... < \j(h) of eigenvalues of 
finite multiplicities, if e and h are small enough. A fortiori, (Hi) insures that E e c 
is compact. 

To simplify notations we write z = (x,£) for any point of the phase space. 

(H 2 ) On E e c , Po has a unique critical point zq = (xo,fo) and near zq : 

N 

Po(z) =E C + Y, pj(z) + 0(\\(z - zo) ir +1 ), k > 2, 

j=k 

where the functions pj are homogeneous of degree j in z — zq. 
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Remark 2 Strictly speaking, one could consider k = 2 with supp(y>) small. But 
there is nothing new here since this case is precisely treated in |3j. 

(. H 3 ) We have (p £ Cq°(R). 

Since we are interested in the contribution of the fixed point z o, to understand the 
new phenomenon it suffices to study a local problem : 

1 r . tE r 

j Zo (E c ,(p,h) = — Tr / e q '~^(f(t)'i/j w (x,hD x )exp(--tP h )Q(P h )dt. (4) 
R 

Here 0 is a function of localization near the critical energy surface E e c and ip £ 
Cq° (T*R”) has an appropriate support near zq. Rigorous justifications are given in 
section 3 for the introduction of Q(Ph) and in section 4 for ip w (x 7 hD x ). The case 
of a critical point which is not an extremum is quite difficult, in particular because 
the singularity is transferred on the blow up of the critical point. To obtain a 
reasonable problem we consider the following hypothesis inspired by Hormander’s 
real principal condition for distributions : 

(Hi) We have Vp*, /0 on the set C(pk) = {9 £ S 2n_1 / pk(9) = 0}. 

Remark 3 Contrary to the case of a local extremum JjjjJ, zq is not isolated on 
E_e c . This imposes to study the classical dynamic in a neighborhood of zq- 

By homogeneity we have Vpfc / 0 on the cone {z £ T*R n \{0} / pk(z) = 0}, 
but mainly (H 4 ) will be used on the unit sphere. C Vk is a smooth manifold of 
codimension 1 which can be equipped with an invariant (Liouville) measure. 

Construction of a geometrical measure. To state the main result clearly, we 
explain how to construct a relevant measure on the line. By ( H 4 ) we can locally 
construct on C(pfc) the 2n — 2 dimensional Liouville form dL via : 

dL( 6 ) A dpk(9) = d0, V 6 » £ C(p k ). (5) 

Here the differential are w.r.t. coordinates on § 2 " -1 , d 6 is the standard surface 
density and we note again pk(0) for the restriction of p^ on § 2n_1 . The form dL 
induces a density on C(pfc) and by continuity we can extend the construction to 
close surfaces p k(9) = e for e small enough. Accordingly we define the integrated 
density as : 

Lvol(ti) = J \dL( 0 )\. ( 6 ) 

{p fc (e)=ii} 

An alternative definition is to compute the volume in § 2n_1 of the pullback : 

= {9 G S ' 1 " 1 / Pk(9) G [0, x}}, (7) 
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and to interpret the result as a measure : 


= J Lvol (s)ds. (8) 

o 

This relation is known in geometry as the co-area formula. The most important 
point, that will be exploited in this work is that (H 4 ) insures that Lvol(u) is smooth 
near the origin. In other words, viewing Lvol(u) as a distribution we obtain that 
0 ^ singsupp(Lvol(u)). 

With these objects the new contributions to the trace formula are given by : 

Theorem 4 Under hypotheses (H 1 ) to ( H 4 ) we obtain the existence of a full 
asymptotic expansion : 


OO 

lz 0 (E c ,ip,h ) ~ h^~ n E E hi log(ft.) m Ay m (y>), (9) 

m— 0,1 j —0 


where the logarithms only occur when ( 2 n + j)/k € N* and A j m G 
As concerns the leading term we obtain : 

(1) If k > 2 n (non-integrable singularity on T,e c ) we have : 

7 Z0 (E c ,(p,h) ~ /iT~ n A 0 ,o(v) + 0(h ?!! £ 1 - n log(h)), as h -»• 0, 
where the distributional coefficient Ao,o( ( / ? ) given by : 

(2n) n k ^( t +~ 1,ip ^) J \Pk{0)\~^ dO + ¥>(*)) J |p fc (0)| ~^d9). 

{Pfc>0} {pfc<0} 

( 10 ) 

(2) If the ratio 2 n/k = q is an integer we obtain logarithmic contributions : 

7z 0 (E c ,<p,h) ~ h^~ n log(h)A 0 ,i(<p) + 0(h^~ n ), as h -> 0 , 

where : 

|irV(i) ' (11) 
R 

(3) For 2 n > k and 2 n/k N the asymptotic is as in 1) with the modified distri¬ 

butions : 


Zn _ -1 

t+ , Fit) 


2n— 


, . Zrx _ 1 

Lvol ) + (tf ,ip(t) 


2n-4? 


( du 2n j ' ' \ du 2n 

where the derivatives w.r.t. w are normalized distributional derivatives. 


Lvol ) , 
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(12) 


The meaning of normalized derivative it that one choose the normalization : 

d 2n _yr 1 d 2n 
dw 2n l\ j ~ x du,2n ’ 
so that the distributional derivatives satisfy : 

^^w 2 ^J(w^=(w^,f(w)), ( 13 ) 

for all / G Cq° with / = 0 in a neighborhood of the origin. With the method we 
employ here this normalization appears naturally in the expansion. The distribu¬ 
tional bracket involving Lvol is detailed in section 6 . 

Remark 5 In cases (1)&(3) the remainders 0(h i T ti ~"log(h)) can be replaced 
by 0(h~ 2 r~~ n ) under the only condition that (2n+l)/k ^ N. An interesting point 
is that the degree k of the singularity on T*M. n has the inverse scaling property on 
the spherical blow-up since the new degree is —2 n/k. 

Results (3)h(2) for q > 2 are not intuitive and are certainly difficult to be reached 
without geometry. In particular one has to work in the dual since both Fourier 
transforms w.r.t. t are distributional. In (3), the order 2 n is arbitrary and the 
result is the same for any derivative of order greater than E(2n/fc). Viewing the 
top order coefficients of the trace as distributions, i.e. : 

lz 0 (E C) <p,h) ~ f(h) ( 7 , 99 ), h -> 0 , (14) 

in all cases at hand we obtain : 

Corollary 6 Under the previous assumptions, singsupp( 7 ) = {0}. 

A similar result presumably holds for all terms of the expansion since the asymp¬ 
totic involves distributions \t\ a log(|t[). Note that Corollary 0 is not obvious in 
view of Eq. 0 - Also, it must be pointed out that results (l)Hz(2) for q = 1 are 
bigger than the standard estimate for non-critical energies for which one obtain : 

h 1 ~ n 

7 (E, <p, h) ~ —— £(0)Vo 1(£ e ), (15) 

where Vol(E s ) is the usual Liouville volume of the regular (compact) surface of 
energy E. Hence the presence of non-integrable singularities on the energy surface 
has a significative spectral effect which can perhaps be measured by eigenfunctions 
estimates as in 0. On the other side, for an integrable singularity, we obtain the 
global result : 

Corollary 7 Under the conditions of Theoremif k < 2n we have : 

h 1 ~ n 

l{E c , (p, h) ~ ^-^■^(0)Vo 1(E £ ; c ), as h 0. (16) 
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This result is a consequence of Theorem Q] and of the results of |3| section 4 to 
which we refer for a detailed proof. Finally, an interesting problem that we have 
not investigated here is the problematic of small h and \E — E c \ estimates for a 
non-integrable singularity onS^ : in this setting Vo1(£_e) of Eg. ifHall diverges as 
E —► E c . See 0 for more details. 

3 Oscillatory representation. 

The construction below is more or less classical. We recall here important facts 
and results concerning the approximation of the propagator by Fourier integral 
operators, or FIO, depending on a parameter h. We recall that : 

7 (E c ,(p,h)= Y y ( Xj ^ Ec ) , 4 = [E c -e,E c + e\, 

*j(h)€l e 

with if G Co°(M) and 1 (/ eo ) compact in T* 1". In this setting the spectrum of 
Ph is discrete in I e for h > 0 small enough and e < £q (see Remark^) and the 
sum is well defined. We localize near E c with a cut-off 0 6 Cq°(]E c — e,E c + e[), 
0 = 1 near E c and 0 < 0 < 1 on 1. The associated decomposition is : 

l{E c , <P, h) = 7i(^ c , <P, h) + 7 2 (E c , f, h), 


with : 

7i (E c ,<p,h)= Y (! E ° ), (17) 

A j{h)el e 

7 2 (E c ,ip,h)= Y Ec )■ ( 18 ) 

A j(h)&I e 

A classical result, see e.g. {5;, is that the sum of Eq. EJ satisfies : 

Lemma 8 71 (E c , tp, h) = 0(h°°), as h —> 0. 

Consequently, asymptotic behaviors of 7 (E c , <p, h) and 72(^01 f, h) are equivalent 
modulo 0{h°°). By inversion of the Fourier transform we obtain : 

Q(P h )f( Ph h E ° ) = ^ J e l ~^0(t)exp(-^tP h )e(P h )dt. 

R 

Since the trace of the left hand-side is exactly 72 (E c , p, h), we have : 

I f ■ tE r i 

7 2 {E c ,ip,h) = —TV J e l ~^~ip(t)ex-p(--tP h )Q(P h )dt. (19) 

R 
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Observe that Eq. m generalizes the Poisson summation formula. 

Let be Uh(t) = exp (—jj-P h ), the evolution operator. For each integer N we can 
approximate Uh(t)@(Ph), modulo 0(h N ), by a Fourier integral-operator, or FIO, 
depending on a parameter h. Let A be the Lagrangian manifold associated to the 
flow of po, i.e. : 

A = {(t,T,x,£,y,r)) G T*R x T*R" x T*R" : r = Po(x,£), (x,£) = $t(y,v)}- 


Theorem 9 Uh(t)Q(Ph) is an h-FIO associated to A. There exist (t) with 

integral kernel in Hormander’s class /(R 2 " +1 , A) and R^\t) bounded, with a L 2 - 
norm uniformly bounded for 0 < h. < 1 and t in a compact subset o/R, such that 
U h (t)Q(P h ) = + h N R[ N \t). 

We refer to Duistermaat 0 for a proof of this theorem. 

Remark 10 By a theorem of Helffer&Robert (Theorem 3.11 and Remark 3.14 of 
d), Q(Ph) is an /i-admissible operator with a symbol supported in p 0 1 (/ e ). This 
allows us to consider only oscillatory-integrals with compact support. 

For the control of the remainder, associated to R^*\t), we use : 

Corollary 11 Let 0i G C“(R) such that 0i = 1 on supp(0) and supp(0i) C I e , 
then \/N G N : 

Tr(Q(P h M^^)) = i-Tr J $(t)ei tE °U™(t)e 1 (P h )dt +0{h N ). 

R 

For a proof of this result, based on the cyclicity of the trace, see [H] or d 

If (xo,£o) ^ ^ and G C°°(R fc x R w ) parameterizes A in a sufficiently small 

neighborhood U of (xo,fo), then for each Uh G I(R fc ,A) and \ G Cq°(T* R fc ), 
su PP(x) C U, there exists a sequence of amplitudes aj = aj(x,d) G Co°(R fc x R N ) 
such that for all L G N : 

X w {x,hD x )u h = VI(a je iv) + 0{h L ). (20) 

—d<j<L 

We will use this remark with the following result of Hormander (see . tome 
4, proposition 25.3.3). Let be (T, t, Xo, £o; Ho, ~Vo) G Afl ow , rj 0 ^ 0, then near this 
point there exists, after perhaps a change of local coordinates in y near y 0 , a 
function S(t,x,r]) such that : 

<t>(t,x,y,r]) = S(t,x,ri) — {y,rf ), ( 21 ) 
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parameterizes Afl ow . In particular this implies that : 

{(*, d t S(t, x , rf),x, d x S(t , x, rj), d v S(t, x, rj), - 77 )} C A flow , 
and that the function S' is a generating function of the flow, i.e. : 


^t{d v S{t,x,rj),rj) = (x,d x S(t,x,rj)). ( 22 ) 

Moreover, S satisfies the Hamilton-Jacobi equation : 

f d t S{t,x,rj) +p 0 (x,d x S(t,x,rj)) = 0 , 

1 S(0,x,O = (x,0 ■ 

Mainly, we will apply this result locally near (xq. £ 0 ) = (yo,Vo), our unique fixed 
point of the flow on the energy surface Y<e c ■ If Co =0 we can replace the operator 
Ph by with £1 7 ^ 0. The spectrum is the same since the new 

operator has symbol p(x, £ — £ 1 , h) with critical point (xo, £ 1 ), £1 7 ^ 0 . 

Hence, for each N £ N* and modulo an error 0(h N ~ d ), the localized trace 
72 (-E’c ; < P, h) of Eq. (Jl9l can be written as : 


7 2 (E c ,<p,h) = J2( 2nh )~ d+J j ei^ s ^ x ^-< x ^ +tE ^ aj (t,x^)ip{t)dtdxd^ 

j<N Ixl 2 ” 

(23) 

To obtain the right power — d of h we apply results of Duistermaat 0 (following 
here Hormander for the FIO, see m tome 4) concerning the order. An h- pseudo¬ 
differential operator obtained by Weyl quantization : 

(2? wh)~% J a(^-^,£)e^ <x_y ’ 4> d£, 

R N 

is of order 0 w.r.t. 1 /h. Since the order of Uh(t)Q(Ph) is — we have : 
if w {x,hD x )U h {t)Q{P h ) ~ ^( 2 nh)- n+ i J a j (t,x,y,r 1 )e^ s ^ x ^~^dy. (24) 


j<N 


R n 


Multiplying Eq. (12411 by <p(t)eh tE ° and passing to the trace we find Eq. 12311 with 
d = n and where we write again aj(t, x, if) for aj(t, x, x, rj). 

To each element Uh of I (R fe , A) we can associate a principal symbol e^ s a pr i nc (uh), 
where S' is a function on A such that £dx = dS on A. In fact, if ip, = I(ae^) then 
we have S = S v = <poi~ l and cr P nnc{uh ) is a section of |A |2 ®M(A), where M( A) 
is the Maslov vector-bundle of A and | A|i the bundle of half-densities on A. If p\ 
is the sub-principal symbol of Ph , in the global coordinates (t,y,rj) on Afl ow the 
half-density of the propagator Uh{t) is given by : 

t 

exp(i J pi($ s (y, -r]))ds)\dtdydr)\i. (25) 

0 
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This expression is related to the resolution of the first transport equation for the 
propagator, for a proof we refer to Duistermaat and Hormander m- 


4 Classical dynamic near the equilibrium. 

The function S of Eq. j2St is related to the classical mechanic and we obtain a 


link with the dynamic generated by po- Precisely, a critical point of the oscillatory 
integral of Eq. (1251) satisfies the equations : 



where the right hand side defines a closed trajectory of the flow inside E e c - Ac¬ 
cordingly, the non-stationary phase method shows that 72(fi c , /i) is asymptot¬ 
ically determined by the closed orbits of the flow on E Ec . We recall that we are 
mainly interested in the contribution of the equilibrium Zo. With ip G Cg° , 

ip = 1 near zq , we write : 



With the additional hypothesis of having a clean flow, the second term can be 
fully treated by the semi-classical trace formula on a regular level. We also observe 
that the contribution of the first term is micro-local and allows to introduce local 
coordinates near Zq . To distinguish the contribution of Zo from eventual closed 
trajectories we use the following result on the dynamic. 

Lemma 12 For all T > 0 there exists a neighborhood Ut of the critical point such 
that 7 ^ 2 for all z G Ut\{zq} and for all t g]-T, 0[U]0, T[. 

For a proof of this Lemma, based mainly on a result of Yorke ESI, see 0. 

With supp(<p) compact we can choose 'ip such that Lemma 1121 holds on s\vpp(jjj) 
for all t G supp(p). Hence, on the support of 'ip there is two contributions : 

1) Points (f, £,£) = (0,x,(;) for (x,£) G E Sc . 

2) Points (t, x,£) = ( t , zq) for t G supp(<p). 

Now, we restrict our attention to the second contribution. In the following, un¬ 
til further notice, the derivatives d are derivatives w.r.t. initial conditions z (resp. 
x, £). If zq is an equilibrium of the vector field X then the linearized flow dexp(t.X)(zo) 
is the flow of the linearized vector field dX(zo). In our setting, the linearized op¬ 
erator is zero and we have : 


d$ t (zo) = exp(O) = Id, Vt. 


(26) 
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Moreover, with ( Hi ) we clearly obtain that : 

<P$t(*o) = 0, Vt, Vj € { 2 ,k - 2 }. (27) 

The non-zero terms of the Taylor expansion of the flow are computed by : 

Lemma 13 Let be Zo an equilibrium of the C°° vector field X and the flow of 
X. Then for all m £ N*, there exists a polynomial map P m , vector valued and of 
degree at most m, such that : 

t 

d m <S> t (z 0 )(z m ) =d$t(zo) I d$- s (z 0 )p m (d$ s (z 0 )( z ),...,d m - 1 <f> s (z 0 )(z m - 1 ))ds. 

0 

See 0 or [B) for a proof. In our setting, the [k — l)-jet of po is flat in z 0 and we 
have Pk-i(yi, —,yk- 2 ) = d fe_1 H po (zo)(yi~ 1 )- H ere i f° r any vector u the notation 
u l stands for (u ,..., u) repeated /-times, the same convention is used below. In view 
of Eq. (I2(il) by integration form 0 to t we obtain : 

t 

d fc - 1 $ t (z 0 )(^" 1 ) = J d fe - 1 H Po (^ 0 )( 2 fe - 1 )c/s = td k - 1 H 9 „{z 0 )[z k - 1 ). (28) 

0 

This provides an explicit formula for the germ of $>t in : 

Mz) = z + + °(IWI fc )- (29) 

We describe now more precisely the singularities of our phase function. Without 
loss of generality, we can assume that zq = 0 

Lemma 14 Near the origin we have : 

S(t,x,£) - (x,£) + tE c = -t(p k (x,£) + R k+1 (x,£) +tG k+ i{t,x,£,)), (30) 

where R k+1 (x,£) = C>(||(a;,£)|| fe+1 ) and G k+1 (t, x,£) = 0(\\(x, £)|| fe+1 ), uniformly 
for t in a compact subset of R. 

Proof. In view of Eq. we search our local generating function as : 

S(t,x, 0 = ~tE c + (x,0 + S k (t,x,0 + 0(||(x,^)|| fc+1 ), 

where S k is homogeneous of degree k w.r.t. Let J be the matrix of the 

usual symplectic form. Comparing terms of degree k — 1 in the implicit relation 
®t{diS(t,x,£),£) = (x, d x S(t, x, £)) provides : 

JVS k {t, x,0 = ^ T y j d fc - 1 $ t (0)((2;,O fc - 1 )- 
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By homogeneity and with Eq. (OHIl we obtain : 

S k (t,x,£) = ^ (( x,^),tJd k ~ 1 H Pk (x,f ) k_1 ) = -tpfc(z,0- 

As concerns the remainders, we have 5(0, x,£) = ( x , £), so that : 

S{t,x,0 - (x,£) = tF(t,x,£), 

where F is smooth in a neighborhood of (x,£) = 0. Now, the Hamilton-Jacobi 
equation imposes that F( 0 ,x,£) = —po{x,f) and we have : 

R k+1 (x,0=po(x,Q -E c - p fc (x,£) = O(||(x,0H fc+1 )- 

Finally, the time dependant remainder can be written : 

S(t, x , £) - 5(0, x, 0 - td t S( 0, x, 0 = 0(t 2 ), 

since by construction this term is of order O (||(x, £)|| fe+1 ) we get the desired result 
when t is in a compact subset of R. ■ 

Remark 15 With Lemma C3 one can compute explicitly terms of higher degree 
for S and <!> t . But we do not need them for the present contribution because of 
some considerations of homogeneity below. 

5 Normal forms and oscillatory integrals. 

Retaining only the coefficient of highest degree w.r.t. h in Eq. (E3, we have to 
study the asymptotic behavior of oscillatory integrals : 

<>= / (31) 

RxT*l“ 

Here, we have temporary discarded the factor (2nh)~ n to avoid it’s constant rep¬ 
etition in the calculations. Since the contribution we study is local, we can work 
with local coordinates and we identify locally our neighborhood of the critical 
point in T*R" with an open of R 2n . With 2 = (x,£) G R 2 ”, we define : 

t,z) = 'L(t,x,£) = 5(t,x,0 - (x,£) + tE c . (32) 

The next Lemma provides a resolution of singularities for w.r.t. C(p k )- 

Lemma 16 Assume Ph satisfies conditions (Hfi) and (Hfi). For all t in a compact, 
after a blow-up w.r.t. z in a neighborhood of Zq, there exists local coordinates p such 
that : 

z) ~ —popi, in all directions where p^ > 0, 

^(f, z) ~ +Pop k , in all directions where pt < 0, 
z) ~ —poViVii locally near C(pfc). 
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Proof. We can assume that zq = 0. To perform the blow-up, we use polar coordi¬ 
nates z = ( r,9 ), 9 G § 2rl-1 (R). By Lemma ITTI near zo we have : 

~ - tr k (p k {6) + rRk+i(r, 9) +trG k +i(t,r,6)), 

~ -tr k (p k (0) +Q{t,r, 9)), 

where p^(0) is the restriction of p^, on S 2ra_1 and Q(t,0, 9) = 0. If pfc(0 o ) ^ 0 and 
t in a compact, then for r < ro we have pfc(0) + Q(t , r, 9) ± 0. We define : 

(r]o,r] 2 ,...,ri 2 n)(t,r, 6 ) = (t, 6 \, 6 2 n -i), 

Vi (t,r,9) = r\p k (0) + Q(t,r,9)\ 

In these coordinates the phase becomes —iioVi if Pfc(^o) is positive (resp. r) 0 r] k for 
a negative value). Near 9 0 , we have : 

^(t, 0,9) = \p k (9)\i?0, Vt, 

hence, the corresponding Jacobian satisfies \ Jr]\(t,Q,9) = |pfc(0)|* ^ 0. 

Now, let 9 0 G C(pfc). Up to a permutation, we can suppose that dg 1 pk(9 Q ) ^ 0. 
We accordingly choose the new coordinates : 

(Vo,Vi,V3,-,V2n)(t,r,9) = (t,r,e 2 ,:.,9 2 n-i), 
m{t,r,9) =p k (0) + Q{t,r,6), 

which are locally admissible since | Jrj\(t, 0, 9q) = |c?6»ipfc(0o)| ^ 0. Finally, lemma 
follows by compactness of C'(pfc) . ■ 

If necessary, we can shrink the support of if) to obtain the existence of the normal 
forms inside supp(</?) x supp(i/>). We define now a partition of unity associated to 
pk . We pick a family of functions <pj G C 00 ^ 2 ™ -1 ) such that : 

C(p k ) C (Jsupp(^), 51 & = 1 near C(p k ). 

j 3 

We can also choose each supp(</; ? ) small enough so that normal forms of LemmalT?)! 
exist in [0,7~o] x supp((/>j). Clearly, this family can be chosen hnite and we obtain 
a partition of unity with = 1 — <f>j . The support of <fo is not connected and 
we define <%, with pfc( 0 ) > 0 on supp( 0 g"), similarly we define <£>q where pfc < 0 , 
so that (fo = 4>o + 4 >o ■ 

Remark 17 The family ((f)j) depends only on p*,, e.g. we can impose = 0 

for |p fc (0)| > e > 0. This point is useful for the globalization in section 6. 
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Let be A = h 1 . We accordingly split up the integral of Eq. m to obtain : 

J±( A) = J e iX9it ' r ' e) 4>t(e)a{t,r6)r^- 1 dtdrdB 

lxl + xS 2 "" 1 

= J e~' lK±rio ' n ^ ) Af(r]o,r] 1 )driodr] 1 = J A^(±Xr]^,r] 1 ), 

RxR + R + 

respectively for the directions where p kid) > 0 and pk(9) < 0. Here the notation A 
stands for the Fourier transform w.r.t. the first argument and is also used below. 
Similarly, the covering of C(pk) gives : 

Ij{ A)= J e iX ' s ( t,r ’ e ^(j)j(9)a(t, r0)r 2n ~ 1 dtdrdQ 

RxR+xS 2 "- 1 

= J e~ lXrioTl ^ V2 A j (7j 0l 7^ 1 ,ri2)dr] 0 dT] 1 d7j2= J A,-( 17 ^ 172 , Vi, V2)dr]idri2- 

Ext+xM R+xR 

These new amplitudes are obtained by pullback and integration, i.e. : 

A o(Vo,Vi) = J r]*{(^{e) a {t,re)r 2n ^ 1 \Jiq\)dri 2 ..-dr] 2 n, (33) 

A j(Vo,m,V 2 ) = J rf((j) j {e)a{t,re)r 2n ~ 1 \J7^\)dr] 3 ...dri2n- (34) 

With C(pk) compact, our oscillatory integral can be written as a finite sum : 

L 

I (A) = 1+ (A) + /_ (A) + ^ Ij (A), A = h~\ (35) 

3=0 

where each term of the r.h.s. will be treated by elementary methods. 

Remark 18 By pullback of the measure r 2n ~ l dr , we have A± = 0(i^l n ~ 1 ) near 
r]i = 0. Same remark for Aj = 0{r^ n ~ 1 ) near 771 = 0. This point is important, 
since Lemmas mm below involve Dirac-distributions w.r.t. 771 . 

Expansion of the related oscillatory integrals. 

We end this section with two results on asymptotics. In fact, to save a lot of 
computations we will take benefit of the linear term rjo in our normals forms. This 
approach is more economic than the strategy proposed in m ■ But the reader must 
keep in mind that the method of m can be applied in a more general setting. The 
next elementary Lemma can be found in [S] and allows to expand both integrals 

I±( A)- 
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Lemma 19 For any a £ Co°(K. x R + ) we have : 



X 


+oo, 


(36) 


where the distributional coefficients are : 

Cj = \— ] {.T(x + k )(t7o) ® 6o j) (m)), x+ = max(x,0). 

AC J ! 

Remark 20 The same result holds for a phase —r/oVi if we replace terms x + by 
x _ in Tiemma.lT!ll For oscillatory integrals on R 2 one can conclude by splitting the 
domain of integration and 2 applications of Lemma l*HH 


For a £ x K + xR), we define the family of elementary fiber integrals : 

OO 

In,k{ A) = J {J a(Xy^y 2 ,yi,y 2 )dy 2 )yl n ~ 1 dy 1 . (37) 

0 R 


Lemma 21 There exists a sequence of distributions (Dj iP ) such that : 

~ E E Dj iP (a) A i log(A) p , as A —> oo, (38) 

p—0,1 j^N, j>2n 

where the logarithms only occur when ( j/k) is an integer. 

As concerns the leading term, if ( 2 n/k) N* we obtain : 

2 TL 2 TL “I - 1 

In,k( A) = A k d(a) + 0 (A k log(A)), (39) 


with : 

OO OO 

d{a) = C n , k J J t^~ 1 yl n ~— [dy™a(t, 0, y 2 ) + d%*a{-t, 0, -y 2 )) dy 2 dt 
0 0 

But when 2n/k = q G N * 7 we have : 

In,k{ A) = A-^log(A)i J \t\ ( i- 1 d^ 1 a(t,0,0)dt + O(X-^). (40) 

R 

Remark 22 The remainder of Eq. cnn can be optimized to 0(X ) if (2 n + 
1 )/k is not an integer and is optimal otherwise. 
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Proof. By a standard density argument we can assume that the amplitude is of the 
form a(s, yi,y 2 ) = f(s)b(y\ , 2 / 2 ). The justification is that our coefficients below are 
computed by continuous linear functionals, i.e. distributions. We define the Melin 
transforms of / as : 

OO 

M±(z) = j s z ~ 1 f(±s)ds. (41) 

0 

We split-up I n) k as J+ and by separating integrations y 2 > 0 and 7/2 < 0. By 
Melin inversion formula, we accordingly obtain : 

J+( A) = 7 ^- J M + (z)X~ z j ( : yiy%)- z b(y 1 ,y 2 )yl n ~ 1 dy 1 dy 2 dz , (42) 

7 R2_ 

where 7 = c + ?’R and 0 < c < fc _1 . Similarly we have : 

J'-(A) = ^ J M_{z)\~ z J (yiy%)~ z b{y 1 ,-y 2 )yl n ~ 1 dy 1 dy 2 dz. (43) 

7 R2_ 

The existence of a full asymptotic expansion is a direct consequence of : 

Lemma 23 The family of distributions z 1 —> ( 2 / 12/2 )~ z on Cq° (R 2 + ) initially defined 
in the domain 9?(,z) < kis meromorphic on C with poles : Zj^ = j/k, j € N*. 
These poles are of order 2 when Zj t k € N* and of order 1 otherwise. 

Proof. We form the Bcrnstein-Sato polynomial bk attached to our problem : 

T{y 2 yi) 1 ~ z := Q^Qyk^y^~ z = b k my 2 y k i)- z , 

k 

bk{z ) = (1 - z) JJ(j - kz). 

3 =1 

If 9?(z) < k~ l , (k + l)-integrations by parts yield : 

J {yiy 2 )~ z f{yi,y 2 )dyidy 2 = ^ J (y 1 y^) 1 ~ z (Tf){y 1 ,y 2 )dy 1 dy 2 . 

R2 R2. 

Now the integral in the r.h.s. is analytic in 3?(z) < 1 + fc _1 . After p iterations the 
poles, with their orders, can be read off the rational functions : 

<«> 

This gives the result since p can be chosen arbitrary large. ■ 
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Hence, the following functions are meromorphic on C : 

0 ± (-) = J {yiy%)~ z Hyu±y2)dyidy 2 - (45) 

A classical result, see e.g. j], is that M±(c + ix) G <S(M X ) when c ^ —N. If we 
shift the path of integration 7 to the right in our integral representation, Cauchy’s 
residue method provides the asymptotic expansion. For any d > c, outside of the 
poles, we obtain that : 

J X~ z M±(z)g ± (z)dz - J X~ z M±(z)g ± (z)dz = ^ res(A ~ z M±g ± )(z j , k )- 

c+iR d+iR c<Zj,k<d 

Since d is not a pole the second integral can be estimated via : 

| J X- z M ± (z)g ± (z)dz\<C(f,b)X- d = 0(X- d ), (46) 

d-\-iR 

where, for each d, the constant C involves the L 1 -norm of a finite number deriva¬ 
tives of b. This will indeed lead to an asymptotic expansion with precise remain¬ 
ders. Applying this method to J+(X) and J_(A) we obtain the existence of a full 
asymptotic expansion of the form : 

4,fc(A) ~ E E C^x-ilogixy. (47) 

p=0,l j'GN* 


Moreover, by Lemma l 2 dl these logarithms only occur when j/k is integer. 

Computation of the leading term. 

To avoid unnecessary discussions and calculations below, we remark that we can 
commute the polynomial weight of Eqs. mm via : 

= b(*)( 2 / 2 ^r^-\ (48) 

k 

b(z) = (1 — z) JJ(j — kz + 2n — 1). (49) 

3 =1 


By iteration, we obtain that the poles are the rational numbers : 

J ~\~ 2 71 — 1 

Zp,j,k,n — P T 7 j j G [1, ..., p G hi. 


By examination of the integrals w.r.t. y\, no residue contributes before : 

2 n 

Z0 = T- 


(50) 
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This result explains the effect of the dimension n and justifies fully Eq. pqi. To 
reach zq we need E(2 n/k) +1 iterations. But, by analytic continuation, any bigger 
integer is acceptable. For the computation of residuum below we use : 

A ~ z M + (z)<B n (z) l (51) 

R 2 

2n—1 1 

*»(--) = n ^zry < 52 > 

1=0 y ’ 


This choice of 2 n iterations is arbitrary but avoids a lot of calculations. 
a) Case of Zo simple pole. 

Here the residue can be computed by the limit z z q. We find : 


lim 




2 77 r 

C ~)M + ( -) J fefo) 


2n -^y 2 1 n - 1 T 2n b(y 1: y 2 )dy 1 dy 2 . 


In particular, we can compute the integral w.r.t. y\ via : 

OO 

J yT^df-^bfy^y^dy, = (2kn - l)!c>^6(0, y 2 ). 
0 


A similar result holds for J_ (A) and we obtain : 


J + (A) = A-TC„ |k M + (y) Jy 2 2 n - 2 -?(d%b)(0,y 2 )dy 2 + R 1 (\), (53) 

0 

oo 

J_(A) = A-TC n , fe M_(^) J yf~^ {d 2 y n 2 b){ 0, -y 2 )dy 2 + R 2 ( A). (54) 

o 

Here C Ut k is the canonical constant : 

1 2n i 

c ^ = lXl—IK' ( 55 ) 

3 =1 J k 

Finally, according to Lemma Ebil each remainder Rj is of order 0{ A - ^ - ) if (2n + 
l)/k N and 0( A -22 ^ 1 log(A)) otherwise. This justifies Remark 12^1 

b) Case of zo double pole. 

If h is meromorphic with a pole of order 2 in £o we have : 

1 d 

res(h)(z 0 ) = - lim —{z - z 0 ) 2 h{z). 

2 Z—>Z 0 oz 
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Applying Leibnitz’s rule to d z (X z M + (z)g + (z)), we obtain : 

J+(A) = BA-Tio g (A) + 0(A-T) ; ( 56 ) 

where the distribution B is computed almost as before. We find : 

OO 

B = ~D n>k M + (^) J y 2 2 n -^(d%b)(0,y 2 )dy 2 , 

0 

2 n 

D n k = lim (z - — ) 2< B n (z). 
k 

But q = 2 n/k is an integer and by integration by parts we obtain : 

OO 

J y 2 2 n - q (d%b)(0,y 2 )dy 2 = (- 1 ) 9 + 1 ( 2 n - q )W~'b{ 0 , 0 ). ( 57 ) 

o 

Since a similar result holds for J_(A), we obtain the desired result by gathering 
all the constants and summation. Finally, we can extend our formulas since all 
coefficients in the expansion are of the form : 

(D\f ®b) = ( D{J } {pi b) , d{ 2 € r>'(K). 

By linearity and continuity, the result holds for a general amplitude a. ■ 

Remark 24 The previous method allows to compute all coefficients of the expan¬ 
sion. In particular the expansion also involves logarithmic distributions (associated 
to double poles). We do not detail all these terms since they have, a priori, no in¬ 
variant formulation in the trace formula. 


6 Proof of the main result. 

We have now the desired results concerning the asymptotic behavior of the trace. 
Hence, to prove the main result, it remains to express the top order coefficients of 
the expansion invariantly. Taking Remark II81 into account, to avoid unnecessary 
calculations we define : 


(t2/i) =vT 1 Ao{t,yi), (58) 

Aj{t,y u y 2 ) =yl n ~ 1 A j {t,y 1 ,y 2 ). (59) 

Note that these definitions have no effect on the Fourier transform w.r.t. t. 
Directions where p k (9) 0. 
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By LemmaETlwe obtain that the first non-zero coefficient is obtained for l = 2n—l 
(see Remark El and is given by 


I I / 2 n — k /n i \ a , \ 1 / 2 n — k 2 \ 

k (9 n _ i)l \l^°l+ fc ® 6 o n ~ ) i A o(vo,Vi)) = j:\\vo\ + k ®6 0 ,A+(iio,m)) ■ 

Since by construction : 


( ? 7o>0) = f a{r)o,0)cl>£(6)\pk{0)\ ^ d6, 

S 2n-1 


(60) 


we obtain that the local contribution, associated to supp(^(J"), is : 

1 / 2 n — k \ C 2 n 

j 0 ( 770 ,0)^> J c/)q (0)\pk(0)\~ — d9, 

S 2n-i 

and a similar result holds on supp(</>Q )■ Now, since a(t,0) = <p{t), cf. Eq. (I251l . the 
directions where p/^f?) ^ 0 contribute as : 


/+(A)- 

'i A ^ ( 

(KI+* - j 

( cf+(d)\p k (9)\-^dd, 

(61) 



§2n 

-1 


I-( A)■ 

4^ 


( ^mpkm~^de. 

(62) 


s 2 »-i 


Microlocal contribution of C( p*,). 

According to the analysis above, we will distinguish out the case fc divides 2 n. 

(1) Case of k > 2 n, integrable singularity on the blow-up. 

Here 2 n/k £]0,1[. According to Lemma I2T1 the contribution of Ij( A) is : 

J |t|^- 1 |i/ 2 |"^ (A j (t,0,y 2 ) + A j (-t,0,-y 2 ))dtdy 2 +R(X). 

R 2 + 

Reminding that y 2 (t,0,0) = pfc(0), we obtain : 

J\y 2 \~^A j (t, 0 ,y 2 )dy 2 = a(t,0) J \p k {6)\~^ 

R+ {pfc(6»)>0} 


Via Eq. (|55|> . by summation /(A) is asymptotically equivalent to : 

\ — — ( r r 

J | p k (9)\-^dd + (tj -1 ,^)^ J \p k {6)\-^dd \ . 

\ {Pfc>0} {Pfc<0} 
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Note that none of the coefficients above are equal unless tp or p*, are symmetric. 
(2) Case of q = 2 n/k integer. 

Here the contribution of each Ij{ A) is dominant since we obtain : 

i j {\) = \\-nog{X) J \tr 1 d^ 1 Aj(t,o,o)dt + o(z-^). 

R 

Unless <7=1, there is no way to take the limit directly, and the geometric properties 
are still hidden in the Jacobian. To reach the result we will use the Schwartz kernel 
technic. Clearly, it is enough to evaluate our derivative and to integrate w.r.t. t. 
With s = (si, S 2 ) € R 2 , we write the evaluation as : 

J e ' l{s ’ ( ' vl ’ V2)) ( is 2 ) q ~ 1 Aj(t,y 1 ,y 2 )dy 1 dy 2 ds. 

Here we have used an oscillatory Schwartz kernel for S Vl <S> Syf 1 - This integral 
representation allows to inverse our diffeomorphism to obtain : 

0 ) = J e i ^ r ^ r ’ e ^\is2) p - 1 a(t, r e)(f j (e)drded^ 

If we extend the integrand by 0 for r < 0, the normalized integral w.r.t. (r, Si) 
provides 6 r . By construction 7/2(0, 9) = p/.(#) and we accordingly have : 

dyf 1 Aj(t, 0,0) = aft, 0)^^y J e i ' tphW (iu) q - 1 <t> j (6)d0du. (63) 

RxS "- -1 

This Fourier integral makes sense with §™ _1 compact. We define here a local 
version of the integrated density : 

Jj(w) = J <j)j(0)dL w (6), (64) 

{p fc (6»)=w} 


where dL w is the density induced by the Leray-form dL Pk : d,pk A dL Pk (0) = dO. 
Note that all these objects can be constructed by mean of local coordinates under 
the only condition that supp (cj>j) is small enough near C'(pfc). Moreover, since p/c 
is continuous on each Jj(w) defines a compactly supported distribution, 

smooth near the origin according to {H 4 ). The sum over all the <p ; j gives the 
geometric contribution : 

-i- J e iuw {iu ) p ~ 1 Y, JjMdwdu = d d J‘™\ 0). (65) 

R2 j 

By integration w.r.t. t we obtain the result in general position. 
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Remark 25 The case 2n = k is directly accessible with : 


/(A) = I^X^ LVol (°)/?(*)*+ 0(* _1 )> 

R 


where LVol(O) is usual the Liouville volume of C(pk)- 

(3) k < 2n and simple pole, non-integrable singularity. 

The distributional coefficients for the positive part of each Ij(X) are : 

OO OO 

{ v j,+ i a ) = C n! k J j t^~ 1 yl n ~ nr (d^A j )(t,0,y 2 )dy2dt. 

o o 

With the same oscillatory technic as above, for pfc(6 l ) > 0 we obtain globally : 


(T+,a) = C n ,k 



d 2n Lvo\{w) 

dw 2n 


J \t\^ l (p{t)dt. 

o 


The duality bracket is well defined since Lvol is compactly supported. A similar 
result holds for the negative part of each Ij(X'). The value of C n ,k, see Eq. (l55l) . 
gives the result stated in Theorem0]by normalizing the distributional derivatives. 
Note that our choice is convenient since d is symmetric. I 

Now we detail the construction of the distributional bracket of the part (3) of 
Theorem QJ Clearly, Lvol is supported in [ inf p fc , sup p&]. Let be y G Cg 0 , 0 < 

S 2 " -1 g2n-l 

X < 1 on R chosen such that y = 1 near the origin and x( u ) = 0 f° r M — £ ; with 
e > 0 small enough. We write the geometric contribution as : 


([/, Lvol) = (U, yLvol) + (U, (1 - x)Lvol). 


Away from the origin, e.g. for u > 0, U is smooth and we obtain directly : 

/ d n 0 „_ 2 n \ I r 2 „ 

Cn,k \du?n U+ “ ’^( w ) Lvol ( u )/ = J. J X(fk{0))\f k (0)\ k d0. 

0 } 

Here the value of C' ra ,/c from Eq. 1551) justifies the normalization of Eq. 1 1 .‘if) . For 

_ 2n 

the singular part of ri + T , we use the local regularity of Lvol(u) near u = 0 and 
integrations by parts to conclude. Finally, using the previous trick we obtain a 
formulation which does not depend on the partition of unity by a double covering 
of the sphere and e small enough. 

Remark 26 The key point here is that we can put in duality the distributions 
d 2n \u\ 2n ~ a and Lvol(«) since their singular supports are disjoints. 
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Effect of the sub-principal symbol. 

Until now we have only considered the case of an operator given by quantization 
of a symbol po- But in presence of a sub-principal symbol pi the construction is 
the same. The only important change, see Eq. il!?51) . is that the amplitude has to 
be modified by : 

t 

a(t,z 0 ) = <p(t) exp(z Jpi($ a (zo))ds). 
o 

But Zq in a fixed point of <b t and a fortiori : 

a(t, 2 0 ) = ip(t) exp(itpi ( 20 )). 

Hence if Pi(zq) = 0, which is the case in many practical situations, the trace for¬ 
mula remains the same. Otherwise, by Fourier inversion formula, the effect is a 
shift on <p by pi(zo) in all integral formulae. 
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metries & Universality in Mesoscopic Systems and IHP-Network, Analysis & 
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